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ABSTRACT 

Planck data have been used to provide stringent new constraints on cosmic strings and other defects. We describe forecasts of the CMB power 
spectrum induced by cosmic strings, calculating these from network models and simulations using line-of-sight Boltzmann solvers. We have 
studied Nambu-Goto cosmic strings, as well as field theory strings for which radiative effects are important, thus spanning the range of theoretical 
uncertainty in the underlying strings models. We have added the angular power spectrum from strings to that for a simple adiabatic model, with 
the extra fraction defined as f\ at multipole i - 10. This parameter has been added to the standard six parameter fit using C0SM0MC with flat 
priors. For the Nambu-Goto string model, we have obtained a constraint on the string tension of Gfi/c 2 < 1.5 x 10~ 7 and f < 0.015 at 95% 
confidence that can be improved to Gji/c 2 < 1.3 x 10" 7 and f < 0.010 on inclusion of high-^ CMB data. For the abelian-Higgs field theory 
model we find, Gyt/ AH /c 2 < 3.2 x 10~ 7 and f m < 0.028. The marginalized likelihoods for f and in the f -£l h h 2 plane are also presented. We have 
additionally obtained comparable constraints on f 10 for models with semilocal strings and global textures. In terms of the effective defect energy 
scale these are somewhat weaker at G/u/c 2 < 1.1 x 10~ 6 . We have made complementarity searches for the specific non-Gaussian signatures of 
cosmic strings, calibrating with all-sky Planck resolution CMB maps generated from networks of post-recombination strings. We have validated 
our non-Gaussian searches using these simulated maps in a P/awc^-realistic context, estimating sensitivities of up to AG/u/c 2 « 4 X 10~ 7 . We have 
obtained upper limits on the string tension at 95% confidence of Gji/c 2 < 8.8 X 10~ 7 with modal bispectrum estimation and Gji/c 2 < 7.8 X 10~ 7 for 
real space searches with Minkowski functionals. These are conservative upper bounds because only post-recombination string contributions have 
been included in the non-Gaussian analysis. 

Key words. Astroparticle physics - cosmology: cosmic background radiation - cosmology: observations - cosmology: theory - cosmology: early 
Universe 
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1. Introduction 

This paper, one of a set associ ated with the 2013 release o f 
data from the P/a^c^] mission ( Planck Collaboration l||2013 ), 
describes the constraints on cosmic strings, semi-local strings 
and global textures. Such cosmic defects are a generic outcome 
of symmetry-breaking phase transitions in the early Universe 
(Kibble 1976) and further motivation came from a potential 
role in large-scale structure formation ( |Zeldovich|l 980; Vilenkin 
1981a). Cosmic strings appear in a variety of supersymmetric 
and other grand unified theories, forming at the end of infla- 
tion (see, for example, Jeannero t et al.|[20Q3] ). However, fur- 



ther interest in cosmic (super-)strings has been motivated by 
their emergence in higher-dimensional theories for the origin of 
our Universe, such as brane inflation. These superstring variants 
come in a number of D— and F- string forms, creating hybrid 
networks with more complex dynamics (see, e.g., |Polchinski| 
2005 ). Cosmic strings can have an enormous energy per unit 
length \i that can give rise to a number of observable effects, in- 
cluding gravitational lensing and a background of gravitational 
waves. Here, we shall concentrate on the impact of strings on the 
cosmic microwave background (CMB), which includes the gen- 
eration of line-like discontinuities in temperature. Comparable 
effects can also be caused by other types of cosmic defects, no- 
tably semi-local strings and global textures. As well as influenc- 
ing the CMB power spectrum, each type of topological defect 
should have a counterpart non-Gaussian signature giving us the 
ability distinguish between different defects, alternative scenar- 
ios, or systematic effects. The discovery of any of these objects 
would profoundly influence our understanding of fundamen- 
tal physics, identifying GUT-scale symmetry breaking patterns, 
perhaps even providing direct evidence for extra dimensions. 
Conversely, the absence of these objects will tightly constrain 
symmetry breaking schemes, again providing guidance for high 
energy theory. For a gener al introduction to cosmic strings and 
other defects, r efer to [V ilenkin & Shellard (2000); Hindmar sh &| 
Kibble| ( [T995] ); |Copeland & Kibble| | |2010) . 

High resolution numerical simulations of cosmic strings us- 
ing the Nambu-Goto action indicate that cosmological networks 
tend towards a scale-invariant solution with typically tens of long 
strings stretching across each horizon volume. These strings 
continuously source gravitational perturbations on sub-horizon 
scales, the magnitude of that are determined by the dimension- 
less parameter: 



~r2 



- 

Impi 



(1) 



where 7] is the energy scale of the string-forming phase transition 
and mpi = ^hc/G is the Planck mass. String effects on the CMB 
power spectrum have been estimated using a phenomenological 
string model and, with WMAP and SDSS data, these estimates 



yield a 2cr upper bound of Gji/c 2 < 2.6 x 10 7 (Battye & Moss 



|2010| ). A consequence is that strings can be responsible for no 
more than 4.4% of the CMB anisotropy signal at multipole i - 
10. 

As we shall discuss, the evolution of Nambu-Goto string 
networks is computationally challenging and quantitative uncer- 
tainties remain, notably in characterizing the string small-scale 



1 Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



structure and loop production. An alternative approach has been 
to use field theory simulations of cosmological vortex- strings. 
These yield a significantly lower number of strings per hori- 
zon volume (less than half), reflecting the importance of radia- 
tive effects on the microphysical scales being probed numeri- 
cally. The degree of convergence with Nambu-Goto string sim- 
ulations is difficult to determine computationally at present, but 
there are also global strings for which radiative effects of com- 
parable magnitude are expected to remain important on cos- 
mological scales. It is prudent in this paper, therefore, to con- 
strain both varieties of strings, labelling the field theory con- 
straints as AH from the abelian-Higgs (local U(l)) model used 
to describe them. Given these quantitative differences, such as 
the lower density, field theory strings produce a weake r con- 
straint G/i A u/c 2 < 5.7 x 10~ 7 using WMAP data alone (|Bevis| 



et al.|2 008 ) and ( [Urrestilla et al.|2011| ). The shape of the string 



induced power spectrum also has a different shape, which allows 
up to a 9.5% contribution at I = 10. These WMAP constraints 
can be improved by adding small-scale CMB anisotropy in a 
joint analysis. The Nambu-Goto strings li mit improves to be- 
come Gfi/c 2 < 1.7 x 10 -7 (using SPT data, |Dvorkin et al.|201l) 
and field theory strings yield G/jah/c 2 < 4.2 x 10" 1 (Urrestilla 
|etal. 12011} . Power-spectrum based constraints on global textures 
were studied in |Bevis et aL] ( |20Q4j ) and ( [Urrestilla eraT]|2008b , 
with the latter paper giving a 95% limit of Gfi/c z < 4.5 x 10" . 



Urrestilla et al. (2008 ) also provide constraints on semi-local 
strings, G///c z < 5.3 x 10" 6 . 

Constraints on cosmic strings from non-Gaussianity require 
high resolution realisations of string-induced CMB maps that 
are extremely challenging to produce. Low resolution small- 
angle and full-sky CMB maps calculated with the full recom- 
bination physics included, have ind icated some evidence for a 
significant kurtosis from strings (Landriau & Shellard 2011). 
More progress has been made increating high resolution maps 
from string lensing after recombination (seejRingeval & Bouchet 
(2012) and references therein) and identifying, in principle, the 
bispectrum and trispectrum, which can be predicted for strings 
analy tically pindmarsh eFaLl[2009l |20T0l |Regan & Shellard 
2010). The first WMAP constraint on cosmic strings using the 
analytic CM B trispectrum yielded G jd/c 2 < 1.1 x 10" 6 at 95% 
confidence ( [Fergusson et al.|2010b| ). An alternative approach is 
to fit pixel-space templates to a map, this meth od was applied to 
global textures templates in |Cruz et al.| ([2007 ) and |Feeney et al.| 
( 2012a|b| ). 



The most stringent constraints that are claimed for the string 
tension arise from predicted backgrounds of gravitational waves 
that are created by decaying loops ( Vilenkin 1981b). However, 
these constraints are strongly dependent on uncertain string 
physics, most notably the network loop production scale and the 
nature of string radiation from cusps, i.e., points on the strings 
approaching the speed of light c. The most optimistic con- 
straint base d on the European Pulsar Timing Array is Gji/c 2 < 



4.0 x 10 9 ( van Haasteren et al.|201 1[ ), but a much mor e conser- 



Sanidas 



vative estim ate of Gfi/c 2 < 5.3 x 10~ 7 can be found in__ 
et al.| ( [20T2] ), together with a string parameter constraint survey 
and an extensive discussion of these uncertainties. Such gravita- 
tional wave limits do not apply to global strings or to strings for 
which other radiative channels are available. 

Alternative topological defects scenarios also have strong 
motivations and we report limits on textures and global 
monopoles in this paper as well. Of particular recent interest are 
hybrid networks of cosmic strings where the correlation length 
is reduced by having several interacting varieties (e.g., F- and 
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D-strings) or a small reconnection probability, p < 1. We expect 
to investigate these models using the Planck full mission data. 

The outline of this paper is as follows. In Sect. 2 we briefly 
describe the different types of topological defects that we con- 
sider, and their impact on the CMB anisotropics . We also dis- 
cuss how the CMB power spectrum is computed and how we 
obtain CMB maps with a cosmic string contribution. In Sect. 3 
we present the defect constraints from the CMB power spectrum 
(with numbers given in Table[2]), while Sect. 4 discusses searches 
for topological defects with the help of their non-Gaussian sig- 
nature. We finally present the overall conclusions in Sect. 5. 



2. Theoretical Modelling and Forecasting 

2.1. Cosmic strings and their cosmological consequences 
2.1 .1 . String network evolution 

A detailed quantitative understanding of the cosmological evo- 
lution of string networks is an essential pre-requisite for mak- 
ing accurate predictions about the cosmological consequences 
of strings. Fortunately, all string network simulations to date 
have demonstrated convincingly that the large-scale properties 
of strings approach a self-similar scale-invariant regime soon af- 
ter formation. If we treat the string as a one-dimensional object, 
then it sweeps out a two-dimensional worldsheet in spacetime 



= AO, 



0, 1, 



(2) 



where the worldsheet parameters £° and f 1 are time-like and 
space-like respectively. The Nambu-Goto action that governs 
string motion then becomes 



(3) 



where y a b = g^ v d a x^diyX v is the two-dimensional worldsheet 
metric (y = det(y^)) induced by the spacetime metric g^. The 
Nambu-Goto action Eq. ^ can be derived systematically from 
a field theory action, such as that for the abelian-Higgs model 
describing U(l) vortex- strings: 



(D tl( />nir<p)-—F, v 



> 2 - 



(4) 

where is a complex scalar field, F^ v is the U(l) field strength 



and = 



+ iAp is the gauge-covariant derivative with e and 



A dimensionless coupling constants. The transverse degrees of 
freedom in cp can be integrated out provided the string is not 
strongly curved, that is, the string width S « tic In <*c L where 
L is the typical radius of curvature. For a cosmological string 
network today with Gyu/c 2 ~ 10" 7 , these two lengthscales are 
separated by over 40 orders of magnitude, so this should be a 
valid approximation. 

In an expanding universe, the Nambu-Goto action Eq. ^ 
yields a Hubble-damped wave equation governing the string 
motion. These equations can be solved numerically, pro- 
vided "kinks" or velocity discontinuities are treated carefully. 
However, they can also be averaged analytically to describe the 
scale-invariant evolution of the whole string network in terms of 
two quantities, the energy density p and the r.m.s. velocity v. Any 
string network divides fairly neatly into two distinct populations 
of long (or "infinite") strings stretching beyond the Hubble 
radius and th e small loops p\ with length / <?c H~ l that the long 
strings create Kibble (1985]). Assuming the long strings form a 



Table 1. Summary of numerical simulation results for the string 
density parameter f defined in Eq. (5]). The Nambu-Goto string 
simu lations a re respectively labelled as MS ( [Martins & She llard 

[20Q6 | 
et al 



6b, RS B (Rin geval et aL] |2007), and BOS ( [Blanco-Pillado 
l.|2011| ). This is contrasted with the much lower density re- 
sults from lattice field theory simulations of vorte x- strings la- 
belled as MMS ( [Moore et al-1|2002> and BHKU JBevis et al. 
|2007b» . 



Epoch 


... MS 


RSB 


BOS 


MSM 


BHKU 


Radiation 


11.5 


9.5 


11.0 


5.0 


3.8 


Matter 


3.0 


3.2 


3.7 


1.5 


1.3 



Brownian random walk characterised by a correlation length L, 
we have 



f2 



and the averaged equations of motion become simply 

2— = 2HUX +v 2 ) + cv, 
dt 

dvoo 

dt 



(5) 



= (l-v2) 



-2Hv 



(6) 



where c measures the network loop production rate and k(v) is a 
curvature parameter with k « 2 V2(l - V2v). This is the velocity- 
dependent one-scale (VOS) model and, with a single parameter 
c, it provides a good fit to both Nambu and field theory simula- 
tions, notably through the radiation-matter transition ( Martins & 
|Shellard|1996) . 

A general consensus has emerged from the three main simu- 
lation codes describing Nambu-Got o string networks ([Martins &| 
Shellard|2006l|Ringeval et al.|2007[[Branco-Pillado et al.|2011| ). 



These independent codes essentially solve for left- and right- 
moving modes along the string using special techniques to han- 
dle contact discontinuities or kinks, including "shock fronting", 
artificial compression methods and an exact solver for piecewise 
linear strings, respectively. The consistency between simulations 
is shown in Table [T] for the string density parameter f defined in 
Eq. ([5}. Averaging yields the radiation era density £ = 10.7 and 
a matter era value f = 3.3. Note that these asymptotic values and 
the intervening matter-radiation transition can be well-described 
by the VOS model Eq. ^ with c = 0.23. The matter era VOS 
value appears somewhat anomalous from the other two simula- 
tions, but this is obtained from larger simulations in a regime 
where convergence is very slow, so it may more closely reflect 
the true asymptotic value. These simulations have also advanced 
the study of string small-scale structure and the loop dist ribution, 
about which there had been less consensus (see, e.g., |Blanco-| 
Pilla do et al.||2011| ). However, note that CMB anisotropy is far 
less sensitive to this issue compared to constraints from gravita- 
tional waves. 

Field theory simulations using lattice gauge techniques have 
also been employed to study the evolution of string networks in 
an expanding universe. Comparatively, these three-dimensional 
simulations are constrained to a lower dynamic range and the 
simulations require the solution of modified field equations to 
prevent the string core width shrinking below the lattice reso- 
lution. On the other hand, field theory simulations include field 
radiation and therefore provide a more complete account of the 
string physics. In Table [T] the lower string densities obtained 
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Fig. 1. The spacetime around a cosmic string is conical, as if 
a narrow wedge were removed from a flat sheet and the edges 
identified. For this reason cosmic strings can create double im- 
ages of distant objects. Strings moving across the line of sight 
will cause line-like discontinuities in the CMB radiation. 




142- 



■ 199 



from two sets of abelian-Higg s simulations are given ( [Moore 
|et al.|2002t[Bevis et al.|2007b) . The evolution can be fitted with 
a VOS model Eq. (|6j) with c = 0.57, which is 150% higher than 
for Nambu-Goto strings. Field theory simulations have further 
important applications, particularly for describing delocalised 
topological defects such as textures, for describing models that 
do not form stable defects like semilocal strings, and because 
they include radiative effects naturally. Radiative effects ob- 
served in current abelian-Higgs simulations are comparable to 
the radiative damping anticipated for cosmological global strings 
and so the AH analysis below should offer some insight into this 
case. 



2.1 .2. String gravity and the CMB 

Despite the enormous energy per unit length yu, the spacetime 
around a straight cosmic string is locally flat. The string has an 
equation of state p z = -p, p x = p y = (for one lying along the 
z-direction), so there is no source term in the relativistic version 
of the Poisson equation V 2 ® = AnG(p+ p x +P y + Pz)- The straight 
string exhibits no analogue of the Newtonian pull of gravity on 
any surrounding matter. But this does not mean the string has no 
gravitational impact; on the contrary, a moving string has dra- 
matic effects on nearby matter or propagating CMB photons. 

The spacetime metric about a straight static string takes the 
simple form, 



Fig. 2. Characteristic CMB temperature discontinuity created by 
a cosmic string. Here, the simulated Nambu-Goto string has pro- 
duced a cusp, a small region on the string that approaches the 
speed of light, which has generated a localised CMB signal. 



where v s is the transverse velocity of the string and y s = (1 - 
v 2 )~ 1/2 . This rather simple picture, however, is complicated in an 
expanding universe with a wiggly string network and relativis- 
tic matter and radiation components. The energy-momentum 
tensor T^x, t) essentially acts as a source term for the metric 
fluctuations that perturb the CMB photons and create tempera- 
ture anisotropies. Essentially, the problem can be recast using 
Green's (or transfer) functions G^ v that project forward the con- 
tributions of strings from early times to today: 

AT C 

— (n, x obs , t ) = I d 4 xG^ y (h, x, x obs , t, t ) 7>(x, t) , 



(9) 



ds 2 = dt - dz 2 - dr 2 - r 2 d6 2 , 



(7) 



where n is the line-of-sight direction for photon propagation 
and x b s is the observer position. The actual quantitative solu- 
tion of this problem entails a sophisticated formalism to solve 
the Boltzmann equation and then to follow photon propaga- 
tion along the observer's line-of-sight. An example of the line- 
like discontinuity signal created by a cosmic string in the CMB 
is shown in Fig. [2] In this case, a string cusp has formed on 
the string, causing a strongly localised signal and reflecting the 
Lorentz boost factor in Eq. (ISl). 



that looks like Minkowski space in cylindrical coordinates, but 
for the fact that the azimuthal coordinate 6 has a restricted range 
< 6 < 2n(l - AGp). The spacetime is actually conical with a 
global deficit angle A = 87rGyU, that is, an angular wedge of width 
A is removed from the space and the remaining edges identified 
(see Fig. [I]). This means that distant galaxies on the opposite 
side of a cosmic string can be gravitationally lensed to produce 
characteristic double images. 

Cosmic strings create line-like discontinuities in the CMB 
signal. As the string moves across the line of sight, the CMB 
photons are boosted towards the observer, causing a relative 
CMB temperature shift a cross the string, given by ( [GottIII|1 985 ; 
|Kaiser&Stebbins|1984| ) 

Y = SnGfiv s7s , (8) 



2.2. Semi- local strings 

The tight constraints on the presence of cosmic strings that we 
will discuss later in this paper start to put pressure on the wide 
class of inflation models that generate such defects ( Hindmarsh 
2011). The power of these constraints would be reduced if the 
strings could be made unstable. This is the basic motivation 
behind semilocal strings: a duplication of the complex scalar 
field cp in the abelian-Higg s action (ft]), occurring naturally in a 
range of inflation models ([D asgupta et al. 2004} | Achucarro et"aL 
2006; Dasgup ta et al.|20 07 ), transforms the stable cosmic strings 
into non- topological semilocal strings ( Vachaspati & Achucarro 
|1991| ) as the vacuum manifold becomes S 5 , which is simply- 
connected. The existence and stability of the semilocal strings 
is thus a question of dynamics rather than due to the topology 
of the vacuum manifold. In general we do not expect to form 
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long strings, but rather shorter string segments, as the semilocal 
strings can have ends. The evolution of these segments is very 
complicated and arises directly from the field evolution, so that 
it is only pra cticable to sim ulate these defects with the help of 
field theory (Urrestilla et al.|2008| ). 



2.3. Global defects 

A large alternative class of defects is due to the breaking of a 
global O(N) symmetry (rather than a gauge symmetry as in the 
case of cosmic strings) of a Af-component scalar field 0. The 
energy density of global defects is significantly less localised 
than those that result from gauge symmetry breaking due to the 
absence of the screening by a gauge field, and there are thus 
long-range forces between the defects. The field self-ordering 
is therefore very efficient for all types of defects with N > 2, 
leading to a generic scaling of the defect energy density with the 
background energy density (see e.g. jDurrer et al.|2002| ). For this 
reason global monopoles (N = 3) do not overdose the Universe 
as their local counterparts would. In this paper we study specif- 
ically the case N = 4 called "texture", which can arise natu- 
rally in many multi-field inflation models that involve a non-zero 
vacuum expectation value and symmetry breaking. In this case 
there are no stable topological defects present, but contrary to 
local texture, global texture can have a non-negligible impact 
on the perturbations in the cosmos, with the field self-ordering 
leading to "unwinding events". In spite of their non-topological 
nature, the field evolution is closely related to the one of lower- 
dimensional stable global defects due to the long-range nature 
of the forces. This is similar to the case of the non-topological 
semilocal strings of the previous section, and indeed the semilo- 
cal example can be seen as an intermediate case between cosmic 
strings and global texture: Starting from the semilocal action, we 
can on the one hand revert to the cosmic string action by remov- 
ing one of the complex scalar fields, and on the other hand we 
find the texture action if we remove the gauge field. 



The normalization of global defects is usually given in terms 
of the parameter s = RnGrj 2 lc 2 when using an action like Eq. ^ 
(with a second complex scalar field but without the gauge fields). 
However, for a simpler comparison with the cosmic string results 
we can recast this in terms of G/^/c 2 = si A and quote limits on 
G///c 2 also for the texture model, as in Urrestilla et al. ( 20Q8| . 



2.4. CMB power spectra from cosmic defects 

The CMB power spectrum from topological defects, like strings, 
is more difficult to compute than the equivalent for inflation- 
ary scenarios that predict a spectrum dominated by an adia- 
batic component with a possible, but highly constrained, isocur- 
vature component. In defect-based scenarios the perturbations 
are sourced continuously throughout the history of the Universe, 
as opposed to adiabatic and isocurvature modes that are the re- 
sult of initial conditions. In principle this requires knowledge of 
the source, quantified by the unequal-time correlator (UETC) of 
the defect stress-energy tensor, from the time of defect forma- 
tion near the GUT scale to the present day — a dynamic range of 
about 10 52 — something that will never be possible to simulate. 
Fortunately, we can use the scaling assumption to extrapolate the 
results of simulations with substantially smaller dynamic range. 
This has allowed a qualitative picture to emerge of the charac- 
teristics of the power spectra from defects, though quantitative 
predictions differ. Here, we will focus on spectra calculated in 



two different ways for cosmic strings, as well as spectra from 
semilocal strings and texture models. 

Defect-based power spectra are dominated by different phys- 
ical effects across the range of angular scales, (i) On large an- 
gular scales the spectra are dominated by an integrated Sachs- 
Wolfe (ISW) component due to the strings along the line-of-sight 
between the time of last scattering and the present day. The scal- 
ing assumption implies that this component will be close to scale 
invariant, although in practice it typically has a mildly blue spec- 
trum, (ii) At intermediate scale the dominant contribution comes 
from anisotropies created at the time of last scattering. In con- 
trast to the strong series of acoustic peaks created in adiabatic 
and isocurvature models, defects produce only a broad peak be- 
cause their contributions are not coherent, (iii) At very small an- 
gular scales, the spectra are again dominated by the ISW effect 
because, rather than decaying exponentially due to the effects of 
Silk damping, there is only power-law decay with the exponent 
being a characteristic of the specific type of defect. 

The standard lore is to treat the defect stress energy tensor, 
Ojjy, as being covariantly conserved at first order, which is known 
as the "stiff approximation". In principle, this means that it is 
necessary to measure two independent quantities from the simu- 
lations, or model them. The other two component are then com- 
puted from the conservation equations. In practice things are a 
little more complicated since it is necessary to provide the UETC 



(10) 



where r is the conformal time and k is the wavenumber. Once 
one has the UETC, then there two ways to proceed. The first 
involves creating realisations of the defect stress-energy whose 
power spectra are computed then averaged to give the total 
power spectrum. The other approach involves diagonalization 
of the UETC. During pure matter or radiation domination, the 
scaling property of defect evolution means that quantities are 
measured relative to the horizon scale, so that the UETC is only 
a function of x = kr and x' - kr' . These functions U(x, x') can 
be discretized and then are symmetric matrices that we can di- 
agonalize. The resulting eigenvectors can be inserted as sources 
into a Boltzmann code, and the resultin g Q are then summed up , 
weighted by the eigenvalues ( Penet al. 1997[ Durrer et al. 2002| ). 
Even though the power spectrum resulting from each "eigen- 
source" exhibits a series of acoustic peaks, the summation over 
many such spectra smears them out, as they are not coherent 
(unlike inflationary perturbations). This smearing-out explains 
why defect power spectra generically are smooth, as mentioned 
above. 

There are also several methods to obtain predictions for the 
UETCs of cosmic strings and other topological defects. The first 
approach we will consider for cosmic strings is to use what 
has beco me known as the Unconnected Segment Model (USM; 
|Albrecht et al.|1997|[T999||Pogosian & Vachaspati|1999| ). In its 
simplest form this models the cosmic string energy momentum 
tensor as that of an ensemble of line segments of correlation 
length £dn(t), moving with an r.m.s. velocity (v 2 ) 1/2 , where du(t) 
is the horizon distance. In addition one can take into account the 
effects of string 'wiggles" due to small-scale structure via a co- 
efficient, J3 = yUeff/yU quantifying the ratio of the renormalized 
mass per unit length to the true value. The model parameters 
(v 2 ) 1/2 and J3 are computed from simulations. In our calculations 
we link the USM sources to the line-of-sight Boltzmann solver 
CMBACT (Pogosian & Vachaspati 1999) to create an ensemble 
of realisations from which we find an averaged angular power 
spectrum. 
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NAMBU 
AH mimic 
AH 




Fig. 3. Cosmic string power spectra used in this analy- 
sis: NAMBU (black dashed), AH-mimic (blue dotted) 
and AH (red solid). The spectra have been normalized 
to equal power at I = 10. The spectra are normal- 
ized the observed WMAP1 value at I - 10 and have 
Gjd/c 1 = 1.17 x 10" 6 , 1.89 x 10" 6 and 2.04 x 10" 6 re- 
spectively. Note that the limits discussed in this paper 
mean that the CMB spectra presented here are less than 
3% of the overall power spectrum amplitude and hence 
the differences observed at high £ do not have a much 
effect. 




Fig. 4. Comparison between global texture (black 
dashed) and semilocal (blue dotted) string power spec- 
tra and the AH field theory strings (red solid), normal- 
ized to unity at I = 10. As expected, the SL spectrum 
lies in between the TX and the AH spectra. The AH 
spectrum was recomput ed for the Planck c osmological 
model with sources from |Bevis e t al. (2010), and the SL 
and TX spectra were taken from Urrestilla et al. (2008 ). 



There are two USM-based models that we will use in this 
analysis which we believe span the realistic possibilities — we 
note a more general appr oach marginalizing ove r three string 
parameters is proposed in Foreman et al.| ( |2011| ) (see also re- 
cent work in |Avgoustidis et al.||2012| ). The first USM model, 
which we will refer to as NAMBU, is designed to model 
the observational consequences of simulations of cosmic string 
simulations performed in the Nambu-Goto approximation. In 
these simulations the scaling regime is different in the radiation 
and matter eras, with (£, (v 2 ) l/2 /c, f3) md = (0.13,0.65, 1.9) and 
(£,<v 2 ) 1/2 /c,/?) mat = (0.21,0.60, 1.5) and the extrapolation be- 
tween the tw o is modelled by using the v elocity dependent one- 
scale model (Martins & Shellard 1996 ). In the second, which 
we will refer to as AH-mimic, we attempt to model the field the- 
ory simulations using the Abelian-Higgs model described below, 
with (£, (v 2 ) l/2 /c,/3) = (0.3, 0.5, 1) independent of time. 



The other approach that we will consider is to measure 
the UETC directly from a simulation of cosmic strings in 
the Abelian-Higgs model, which we will refer to as AH. The 
Abelian Higgs model involves a complex scalar field and a 
gauge field described earlier Eq. ([4]), for which the dimen- 
sionless coupling constants e and A are chosen with A = 2e 2 , 
so that the characteristic scales of the magnetic and scalar en- 
ergies are equal, (see |Bevis et al. ( 2007b|a ) for further details 
about the model). We then simulate the evolution of the fields 
on a grid, starting from random initial conditions designed to 
mimic a phase transition, followed by a brief period of diffusive 
evolution, to rapidly reach a scaling solution expected to be typ- 
ical of the configuration found long after the phase transition. As 
the simulation is performed in comoving coordinates, the string 
width is effectively decreasing as time passes. To enlarge the dy- 
namical range available, we partially compensate this shrinking 
with an artificial string fattening. We perform runs for various 
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values of the fattening parameter to ensure that the results are 
not affected by it. 

During the simulations, we compute the energy-momentum 
tensor at regular intervals and decompose it into scalar, vec- 
tor and tensor parts. We store these components once scaling 
is reached, and compute the UETCs by correlating them with 
later values of the energy-momentum tensor. UETCs from sev- 
eral runs are averaged, diagonalized and then f ed into a mo di- 
fied version of the CMBEASY Boltzmann code (Doran 2005]) to 
compute the CMB power spectra (both temperature and polar- 
ization). The spectra used in this paper were derived from field- 
theory simulations on a 1024 3 grid and used the extrapolation 
to sub-string scales described in |Bevis et al.| ([2010), which are 
expected to be accurate at the 10% level to ^ max « 4000. 

In Figs. [3] and [4] we present the spectra we will use in sub- 
sequent analysis. The higher dashed black curve is the spec- 
trum computed using the USM for the NAMBU model, and the 
smaller dashed blue and solid red curves the AH-mimic model 
and the AH model, respectively. We should note that when nor- 
malized to the amplitude of the observed CMB anisotropies 
on large-scales at I = 10, the three models give G/u/c 2 = 
1.17 x 10" 6 , 1.89 x 10" 6 and 1.9 x 10" 6 for the NAMBU, AH- 
mimic and AH models, respectively. The reasons for differences 
between the spectra for these two approaches are discussed in 



Battye & Moss (2010). Briefly, the main reasons for the differ- 



ences are twofold: First, the overall normalization, which is due 
to the NAMBU models having smaller values of more strings 
per horizon volume, and larger values of with each of the 
string segments being heavier, than the two AH models. Both 
these effects mean that a lower value of G/i/c 2 is required to 
achieve the same amplitude for the anisotropies. Secondly, the 
enhanced peak at small angular scales, which is caused by the 
value of £ being smaller in the radiation era than in the matter 
era, meaning that there are more strings per horizon volume in 
the radiation era when the small-scale anisotropy is imprinted, 
and hence more anisotropy on those scales for a given G/u/c 2 . 

The method used for the semilocal strings (denoted SL) and 
0(4) global texture (denoted TX) is fundamentally the same as 
for the AH model: we simulate the field theory on a discretized 
grid and compute the energy-momentum tensor at regular inter- 
vals. From these snapshots we derive the UETCs by correlating 
the scalar, vector and tensor parts at different times. The only 
difference is the field-theory action being used in the simula- 
tions. In Fig. [4] we present the spectra we used for the semilocal 
strings and global textures, taken from Urresti lTa et al.| ( |2008| ). 
These models are also shown with the AH cosmic string model 
for comparison. 

2.5. Maps of CMB anisotropies from cosmic strings 

In order to go further than the two-point correlation function, we 
have used numerical simulations of Nambu-Goto cosmic string 
evolution in an FLRW spacetime to generate various CMB syn- 
thetic maps. The use of simulations is crucial to produce real- 
istic string configurations on our past light cone and have been 
the subject of various code development in the last twenty years 
(see|Albrecht & Turok|ll9"89l |Bennett & Bouchet|[T989l JT990} 
Allen & Shellard 1990; Vincent et al.| 19981 |Moore et al.|200"Tt 
Ringeval "et al.|2007| |Blanco-Pillado et al.|2011| ). Until recently, 
the underlying numerical challenges have limited the resolu- 
tion of the full sky maps to an angular resolution of 1 ^(corre- 
sponding to a HEALPix resolution of N s ^ e = 256) in Landriau 
& Shellard| ( [2003] [20TT] ) (see also early work in |Allen et al. 
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Fig. 5. Integrated Sachs-Wolfe angular power spectra extracted 
from the full sky cosmic string maps at different resolutions (la- 
belled by A^ S ide), with or without applying the anti-aliasing proce- 
dure (see text). The anti-aliasing filtering gives back the correct 
power up to 4iax < 2Af side . 



s mall scales probed b y Planck, we have used the maps described 
in [Ringeval & Bouchet] ( |2012| ) that have an angular resolution 
of 0.85' (A^ S ide = 4096). This map is obtained by consider- 
ing the ISW contribution from (9]), sourced by the Nambu-Goto 
stress tensor, and w hich can be recast into the form ( Stebbi ns~&| 
Veeraraghava nj 1 995] ) 



yW 



4GyU 



f 

Jxn* 



(h • X') • X' 
l+h-X 



Xh-X 
(Xn - X) 2 



dl. 



(11) 

The integral is performed over all string position vectors X = 
{X 1 } intercepting our past line cone (in the transverse temporal 
gauge). Primes and dots denote differentiation with respect to 
the spatial and time-like worldsheet coordinates £ 1 and £° re- 
spectively, while dl is the invariant string length element. Taking 
the limit Xh — > X gives back the sma ll angle and flat sky ap - 

<T588); 



proximation used in |Hindmarsh (1994); Bouchet et al 



Frai sse et al.| ( |2008) . For generating the full sky map, Eq. ([TT 
has been evaluated without any other approximation and re- 
quired more than 3000 Nambu-Goto string simulations to fill the 
whole comoving volume between the observer and the last scat- 
tering surface. De tails on the numerics can be found in Ringeval 
|& Bouchet| ( |20T2] ). 

This method therefore includes all string effects from the last 
scattering surface till today, but does not include the Doppler 
contributions induced by the strings into the plasma prior to 
recombination. As a result, our full sky map represents the 
ISW contribution from strings, which is dominant at large and 
small scales but underestimates the signal on intermediate length 
scales, as can be directly checked on the power spectrum (see 
Fig. [5]). We therefore expect the string searches based on this 
map to be less constraining than those using the power spec- 
trum, but certainly more robust as any line-like gravitating object 
should generate such a signal. 

Calibration and training for the non-Gaussian searches of 
Sect. [4] have required the generation of new full sky and sta- 
tistically independent cosmic string maps. The numerical chal- 
lenges underlying the A^de = 4096 map (Ringeval & Bouchet 
2012) are such that it was numerically too expensive to create 



1996 ). In order to extend the applicability of these maps to the another one of the same kind. At this resolution, the computa 
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Fig. 6. All sky Mollweide projection of the simulated cosmic strings CMB sky after convolution by a Gaussian beam of 5' resolution. 
The color scale indicates the range of AT/T/(G/u/c 2 ) fluctuations. 




-80.0 ^^^H l^^^m 80.0 ^^^m 25000 

AT/T/(G^/c 2 ) |V(AT/T)|/(G / x/c 2 ) 



Fig. 7. A 20° gnomic projection patch extracted from the full sky map and zooming into string induced temperature steps (see 
Fig. [6]). Applying the spherical gradient magnitude operator enhances the temperature steps, and thus the string locations, even 
more (right). 



tions typically require 800000 cpu-hours, so we have chosen 
to generate three new maps at a lower resolution of 1.7', i.e., 
N S ide = 2048. Unfortunately, at this lower resolution, the sim- 
ulated string maps, hereafter referred to as raw maps, exhibit 
a strong aliasing at small scales that could have induced spuri- 
ous systematics even after convolution with the Planck beam. 



This aliasing concerns pixel-sized structures and comes from 
the method used to numerically evaluate Eq. ( pTj ). In order to 
spare computing time, the signal associated with each pixel is 
only computed at the centroid direction h. This has the effect of 
including some extra power associated with string small-scale 
structure that is below the pixel angular size, thereby aliasing 
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the map. In order to address this problem, we have used semi- 
analytical methods to design an optimal anti-aliasing filter, both 
in harmonic space and in re al space. As discussed in |Fraisse etal.| 
(2008); |Bevis et al.| ( |2010| ), the small scale angular temperature 
power spectrum slowly decays as a power law £~ p such that any 
deviations from this behaviour can only come from the aliasing. 
For each N^te = 2048 raw map, we have performed a multi- 
parameter fit of the power spectrum, and of the one-point dis- 
tribution function, to extract, and then removes, its small scale 
aliasing contribution. In order to validate the procedure, we have 
checked that the power spectrum of each of the filtered maps 
matches the one associated of the raw Af^e = 4096 map, the lat- 
ter being also being affected but at half the scale. In Fig. [5] we 
have plotted the power spectra of one of the 7V S id e = 2048 maps 
before and after convolution with our anti-aliasing filter. As ex- 
pected, it matches with the one extracted from the A^e = 4096 
map (here truncated at £ = 4096). Finally, in order to include 
the effects associated with the HEALPix pixelization scheme, the 
anti-aliased maps have been convolved with the HEALPix pixel 
window function before being used for further processing. 

In total, this method has provided four theoretical full sky 
string maps that have been used in the string searches we will 
discuss in Sect. 4] As an illustration, we have represented in 
Fig. [6] one of the filtered string map after convolution by a 
Gaussian beam of FWHM = 5'. The color scale traces the rel- 
ative temperature anisotropies AT/T, divided by the string ten- 
sion G///c 2 . The anisotropy patterns may look Gaussian at first 
because most of the string signatures show up on the smallest 
length scales. In Fig.[7j we have plotted a gnomic projection rep- 
resenting a field of view of 20°, in which the temperature steps 
are now clearly apparent. The right panel of Fig.|7]represents the 
magnitude of the spherical gradient, which enhances the steps. 

Finally, in order to provide a much larger statistical sample 
beyond only four string realisations, we have also produced a 
collection of 1000 small angle patch es (7.2') of the CMB sky de- 
rived in the flat sky approximation (Stebbins 1988; Hindmarsh 



1994} |Stebbins & V eeraraghavan| |1995[ |Bouchet et al.||1988[ 



Frai sse et al.||2008| ). Although the large-scale correlations are 



lost, these maps have been shown to accurately reproduce var- 
ious analytically expected non-Gaussian string effects such as 



the one-point and higher ^-points functions by Takahashi et al. 


(2009};|Hindmarshetal.|( 


2009 2010); Regan &Shellard (2010); 


Yamauchi et al. (2010b|a 


);Ringeval(2010). 



3. Power spectrum constraints on cosmic strings 
and other topological defects 

In order to compute constraints on cosmic string scenarios we 
just add the angular power spectrum to that for an simple adia- 
batic model — which assumes that they are uncorrelated — with 
the fraction of the spectrum contributed by cosmic strings be- 
ing fio at £ = 10. This parameter is then added as an ex- 
tra parameter to the standard six par ameter fit using C0SM0MC 
and the Planck likelihood described in Planck Collaboration XV 
( [2013] ). We use a Flat ACDM cosmology defined through the 
physical densities of baryons, fl^h 2 , and cold dark matter, Q c /z 2 , 
the acoustic scale, #mc, the amplitude, A s and spectral index, 
n s of density fluctuations and the optical depth to reionization 
r. The Hubble constant is a derived parameter and is given by 
Hq = 100 /z km sec" 1 Mpc -1 . We use the same priors on the 



cosmological and nuisance parameters as are used in Planck 
Collaboration XVI| ( |2013| ) and use WMAP polarization data to 



Table 2. 95% upper limits on the constrained parameter fio 
and the derived parameter G/i/c 2 for the five defect models 
discussed in the text. We present limits using Planck and po- 
larization information from WMAP (Planck + WP), and from 
also including high £ CMB information from ACT and SPT 
(Planck +WP+highL). 



Defect type 


Planck+WP 


P/arcc£+WP+highL 




• • • fio 


Gii/c 2 


fio 


Gia/c 2 


NAMBU 


0.015 


1.5 x 10" 7 


0.010 


1.3 x 10" 7 


AH-mimic .... 


0.033 


3.6 x 10" 7 


0.034 


3.7 x 10" 7 


AH 


0.028 


3.2 x 10" 7 


0.024 


3.0 x 10" 7 


SL 


0.043 


ll.Ox 10" 7 


0.041 


10.7 x 10" 7 


TX 


. . . 0.055 


10.6 x 10" 7 


0.054 


10.5 x 10" 7 



added hig h-l CMB data from SPT and ACT t o obtain stronger 
constraint( |Sievers et al.|20131 |Hou et al.|20T2) . 

For the USM-based models we use the approach used in 
|Battye et al.] p006) and |Battye & Moss] pOlt)) . We find that 
the constraints on the standard six parameters are not signif- 
icantly affected by the inclusion of the extra string parameter 
and that there are no significant correlations with other param- 
eters (see Table 13). For the case of Planck data only and us- 



ing the NAMBUmodel we find that Gjd/c 2 < 1.5 x 10" 7 and 
fio < 0.015, whereas for the AH-mimic model we find that 
Gjd/c 2 < 3.6 x 10 -7 and f i0 < 0.033, with all the upper lim- 
its being at 95% confidence level. The ID marginalized likeli- 
hoods for fio are presented in the upper panels of Fig. [5] The 
differences between the upper limits for the NAMBU and AH- 
mimic models is compatible with those seen previously using 
WMAP 7-year and SDSS data ( |Battye & Moss|2010| ). The upper 
limits from this version of the Planck likelihood are better than 
those computed from WMAP7+SPT ( [Dvorkin et luTpOTT) and 
WMAP7+ACT ( [Dunkley et al.|201 1|) and are significantl y better 
than those from WMAP7+SDSS ( [Battye & Moss|20T0) . Based 
on the Planck "Blu e Book" values for noise levels we predicted 
( [Battye etal.|2008) a limit of Gji/c 2 < 6x 10" 8 , while the present 
limit is about a factor of two worse than this. The main reason 
for this is that the projected limit ignored the need for nuisance 
parameters to model high £ foregrounds and that not all the fre- 
quency channels have been used. The corresponding limits for 
the AH model are f i0 < 0.028 and G/i/c 2 < 3.2 x 10" 7 . 

There is now very little degeneracy between the fio and 
parameters, something that was not the case for WMAP alone 
( [Battye et al]|2006l |Bevis et al.|[2QQ8] |Urrestilla et"aLl|20TT) . 
This has implication for supersymmetric hybrid inflation mod- 
els as discussed in Batt ye et al.| ( [201 0\ that typically require 
fts > 0.98. The simplest versions of these models appear to be 
ruled out. The strongest correlation using the NAMBU and AH 



10 



help fix r. In addition to just using the Planck data, we have also 



mimic models is between / 10 and tl^h 2 as illustrated in Fig. 
In addition, we find agreement with Lizarrag a et al. (2012), that 
there are significant correlations between the amount of strings 
fio in the AH model and the number of relativistic degrees of 
freedom 7V e ff as well as between fio and the primordial helium 
abundance Yu e . We leave a detailed study of these correlations 
to later work. 

In Fig. [8] we also present the ID marginalized likelihoods 
for the texture and semilocal string models (compared to the 
AH field theory strings). The resulting constraints on the fio 
parameter are given in Table [2] as well. For the conversion 
into constraints on G/u/c 2 we have that for semilocal strings 
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Table 3. Constraints on the fitted cosmological parameters in the case of Planck alone. It is clear from this that the fitted parameters 
are not significantly affected by the inclusion of defects. 



Parameter NAMBU 

Q h h 2 0.0223 + 0.0003 

Q c h 2 0.119 + 0.003 

<9 MC 1.0415 ±0.0006 

r 0.089 + 0.013 

log(10 10 A s ) 3.080 + 0.027 

n s 0.961 + 0.007 

H 68.4 + 1.3 

Gid/c 2 <1.5xl0" 7 

/io < 0.015 



AH mimic 



AH 



SL 



TX 



0.0223 + 0.0003 
0.119 + 0.003 
1.0415 ±0.0006 
0.090 + 0.013 
3.080 ± 0.026 
0.963 ± 0.008 
68.3 ±1.2 
< 3.6 x 10" 7 
< 0.033 



0.0223 ± 0.0003 
0.119 ±0.003 
1.0415 ±0.0006 
0.090 ± 0.013 
3.081 ±0.025 
0.963 ± 0.008 
68.3 ± 1.3 
< 3.2 x 10" 7 
< 0.028 



0.0223 ± 0.0003 
0.119 ±0.003 
1.0415 ±0.0006 
0.090 ±0.013 
3.081 ±0.025 
0.964 ± 0.007 

68.2 ± 1.2 
< 1.10 x 10" 6 
< 0.043 



0.0223 ± 0.0003 
0.119 ±0.003 
1.0415 ±0.0006 
0.088 ± 0.014 
3.078 ± 0.028 
0.965 ± 0.008 

68.3 ±1.2 
< 1.06 x 10" 6 
< 0.055 
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Fig. 8. Marginalized constraints on fio for topological defects from Planck data plus polarization from WMAP (Planck+WP). The 
left panel show constraints on cosmic strings, with NAMBU in black dashed, AH-mimic in blue dotted and AH in red solid. The 
right panel show the constraints on SL (blue dotted) and TX (black dashed) compared to AH (again solid red). 
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Fig. 9. Marginalized constraints on fio for topological defects with high-/' CMB data from SPT and ACT added to the Planck + 
WP constraints data (compare with constraints shown in Fig.[8|. The left panel show constraints on cosmic strings, with NAMBU 
in black dashed, AH-mimic in blue dotted and AH in red solid. The right panel show the constraints on SL (blue dotted) and TX 
(black dashed) compared to AH (solid red). 



Gfiio/c 2 = 5.3x10 6 and for global texture Gyu io /c 2 = 4.5x10 6 , Gyu, semilocal strings lead to significantly less anisotropies than 
cf |Urrestilla et al.| p008 ). We notice that, as expected for a fixed cosmic strings (a factor of about 8 in the Q), and texture are sim- 
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Fig. 10. Marginalized likelihoods for f\o and in the fio-fl^h 2 
plane for the NAMBU model in blue and the AH mimic model 
in red using Planck +WP. This is the strongest correlation with 
any of the standard cosmological parameters 

ilar to the semilocal strings. We thus expect significantly weaker 
constraints on G\i for the SL and TX models, especially since 
in addition the constraints on /io for these models are weaker. 
Indeed we find a 95% limit of G/i/c 2 < 1.10 x 10" 6 for semilo- 
cal strings and G/u/c 2 < 1.06 x 10" 6 for global textures. 

4. Non-Gaussian searches for cosmic strings 

Cosmic strings and other topological defects generically create 
non-Gaussian signatures in the cosmic microwave sky, counter- 
parts of their inevitable impact on the CMB power spectrum. 
This is a critical test of differentiating defects from simple infla- 
tion, while offering the prospect of direct detection. Searches for 
these non-Gaussian defect signatures are important for two key 
reasons: on the one hand, constraints from the CMB power spec- 
trum can be susceptible to degeneracies with cosmological pa- 
rameters in the standard concordance model; on the other hand, 
any apparent defect detection in the power spectrum should have 
a well-defined prediction in higher-order correlators or other 
non-Gaussian signals, and vice versa. Non-Gaussian tests can 
also be used to distinguish cosmic defects from residual fore- 
grounds or systematic contributions. Below we will present re- 
sults from NG tests that seek strings in multipole space (bis- 
pectrum) and in real space (Minkowski functionals), as well as 
hybrid methods (wavelets). 

4.1. Foregrounds, systematics and validation 

It is well-known that the microwave sky contains not only the 
CMB signal but also emission from different astrophysical con- 
taminants. In particular, point source emission is expected to be 
a special cause of confusion for cosmic defects, notably those 
with high resolution signatures, such as cosmic strings. In addi- 
tion, systematic effects may also be present in the maps at a cer- 
tain level. Therefore, before claiming a cosmological origin of a 
given detection, alternative extrinsic sources should be investi- 
gated and discarded. This can be done by performing a number 
of consistency checks in the data, most of which are common 
to the other non-Gaussianity papers, where they are discussed 
in greater detail. Here, we provide a brief summary of the main 
issues. 



Foreground-cleaned CMB maps are provided using four dif- 
ferent component separation techniques (for further details, see 
Planck Collaboration XII 2013 ): SMICA (semi-blind approach); 
NILC (internal linear combination in needlet space) ;SEVEM (in- 
ternal template fitting); and Commander/Ruler (C-R, parametric 
method). To demonstrate the robustness of a particular result, 
it should be replicated with at least two different cleaned CMB 
maps. The adoption of different masks that exclude different re- 
gions of the sky (ranging from more aggressive to more conser- 
vative) has been used to test the stability of non-Gaussian es- 
timates. Another interesting test is the use of cleaned maps at 
different frequencies (for instance, those provided by the SEVEM 
foreground separation technique). A given detection should be 
consistent at all frequencies, since the behaviour of contaminants 
and systematic effects will, in general, vary with frequency. A 
further test is the study of noise maps constructed from the dif- 
ference between two Planck maps (either at the same or at dif- 
ferent frequencies) smoothed to the same resolution. These maps 
will not contain the CMB signal and, therefore, any NG detec- 
tion should vanish on them. The opposite would indicate that the 
claimed result is due to foreground residuals or to the presence 
of systematic effects. 

In order to validate the various non-Gaussian methods 
that are described below, we instituted a series of Planck 
String Challenges. These were blind tests employing post- 
recombination string simulation maps with an unknown Gyu/c 2 , 
co-added to a Gaussian CMB map, together with standard masks 
and increasingly more realistic noise models. For calibrating the 
non-Gaussian tests, several different string simulations were also 
provided. In addition, 1000 ACDM Gaussian maps, with simu- 
lated noise and using the same mask, were provided for anal- 
ysis purposes, notably for determining the variance of different 
techniques. The aim has been to determine the sensitivity of the 
proposed non-Gaussian tests and to see if the G/i/c 2 in the string 
challenge map could be measured accurately. The results from 
these challenges were an important part of the validation for each 
of the methods described below. Planck simulation pipelines for 
each of the component separation methods were also used to es- 
timate realistic foreground residuals and were used to validate 
non-Gaussian pipelines and to remove string signal bias. 

4.2. Cosmic string bispectrum 
4.2.1. Modal bispectrum methods 

The CMB bispectrum is the three point correlator of the a\ m coef- 
ficients, ^fm 2 m 3 - U{\m x a£ 2 ,n 2 a hnn- F° r the purposes of a search 
for cosmic strings we assume the network cumulatively creates 
a statistically isotropic signal, that is, we can employ the angle- 
averaged reduced bispectrum b^ 2 ^ 3 , defined by 

777/ 

where is a weakly scale-dependent geometrical factor and 
Qmjmlm * s me well-known Gaunt integral over three Y^ m s that 
can be expressed in terms of Wigner-3 j symbols. The CMB bis- 
pectrum b^ 2 f 3 is defined on a tetrahedral domain of multipole 
triples (A satisfying both a triangle condition and i < £ max 
set by the experimental resolution. When seeking the string bis- 
pectrum in the Planck data, we employ the following es- 
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timator to find or limit its amplitude: 



1 0m\ni2iu\ 



string 



(13) 



where we assume a nearly diagonal covariance matrix 
Q imi ,£ 2m2 ~ Ci8 ili2 6 mi - m2 and we modify Q and b hhh ap- 
propriately to incorporate instrument beam and noise effects, as 
well as a cut-sky. To simplify Eq. ( 13 ), we have ignored the "lin- 
ear term" (which is included in the analysis). A much more ex- 
tensive introduction to bispectrum estimation can be found in 
IPlanck Collaboration XXIVl ( |20T3] ). 

A key step in observational searches for non-separable bis- 
pectra, such as those induced by cosmic strings (denoted by 
Z?^ r ), is to expand it into separable modes (Fergusson & 



Shellardl [20091 |Fergusson et al.||2010a"| taking the signal-to- 
noise-weighted form: 



.string 



(14) 



where the modes Q n (tuh, Zi) = \\4 P (h) qSi) q s (h) + perms.] 
are constructed from symmetrized products (the n label a 
distance-ordering for the triples {prs}). The product basis func- 
tions Qnif i^i^y) are not in general orthogonal, so it is very 
useful to construct a related set of orthonormal mode functions 
Rnifiuti* £3) sucri that (R n , R p ) = S np . Substituting the separable 
mode expansion ( [T4] ) reduces the estimator (13) to the simple 
form 



N 2 



(15) 



where the ft® coefficients are found by integrating products of 
three Planck maps filtered using the basis function a% (an effi- 
cient product with each map multiplied by the separable q r (£)). 
We can validate this estimator by using the modal methodology 
to create CMB map realisations fo r cosmic strings from t he pre- 
dicted ft® with a given Gfi/c 2 (seelFergusson et al. |2010a| . It is 
easy to show that the expectation value for J3" for such an en- 
semble of maps in the orthogonal basis should be 



(16) 



Alternatively, we can exploit this fact by reconstructing the a R n 
from given CMB map realisations created directly from string 
simulations, an approach we will adopt here. 

4.2.2. Post-recombination string bispectrum 

In order to estimate the string bispectrum at Planck resolu- 
tion we employed the modal reconstruction method Eq. (16) on 
the post-recombination string challenge simulations described in 
Sect 2.5. These string maps include the accumulated line-like 
discontinuities induced by the string network on CMB photons 
propagating from the surface of last scattering to the present day. 
This work does not include recombination physics, that is, con- 
tributions from the surface of last scattering that will increase the 
string anisotropy signal substantially. As discussed previously, 
there are four full-sky maps of two different resolutions, which 
were provided for the purpose of calibrati ng Planck searches for 
cosmic strings (Ringeval & Bouchet 2012). For the modal anal- 
ysis, we have adopted the hybrid polynomial basis augmented 
with a local shape mode (in total with n max = 600 modes), as 
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Fig. 11. Coefficients a R n Eq. ( 14 ) for the hybrid Fourier mode ex- 
pansion of the cosmic string bispectrum Eq. ( [T4] ). The average 
value (a%) (bold black line) is in remarkable agreement with 
all four string simulations as can be seen for n < 30 (lower 
panel), with each exhibiting better than a 97% correlation over- 
all. Figure to be updated imminently. ] 



well as a hybrid Fourier basis (n mdLX = 300), which are both de- 
scribed in |Planck Collaboration jOOVl ( [20T3] ). 

To extract the string bispectrum in a Planck-realistic context, 
we chose a fairly high non-Gaussian signal with Gji/c 2 = 1 x 
10" 6 . The normalized string maps were added to noise maps gen- 
erated by the component separation pipelines of SMICA, NILC 
and SEVEM, creating twelve sets of 200 simulated string realisa- 
tions. Each of these maps was then filtered using the modal esti- 
mator to find the J3„ coefficients appropriate for each component- 
separation method. After averaging each set of modal coeffi- 
cients a% = (fi^) over the different (unmasked) noise realisa- 
tions, we found remarkable consistency between the estimated 



P„ for the four string simulations as shown in Fig. 11 for the 



Fourier modes. Agreement was good across all the 300 a% 
modes determined, as shown in detail for n=l-30 (see the lower 



panel of Fig. 11 



Quantitatively, the different string simulations produced bis- 
pectrum shapes that had above 97% correlations with each 
other (i.e., Z n N~ l a ( n l) > 0.97 for N 2 = X„oi 1)2 Z^ 2 )- 
The overall integrated bispectrum amplitudes was consistent to 
within 4%. Despite only four string map simulations, these are 
small errors relative to experimental and theoretical uncertain- 
ties. This robustness indicates that the overall string bispectrum 
signal at Planck resolution is a statistical summation of very 
many contributions from the millions of strings between the ob- 
server and the last scattering surface. To ensure the bispectrum 
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C/ weighting was not significantly affected by the presence of 
the large string signal, we repeated the modal extraction pro- 
cedure for Gyu/c 2 = 5 x 10" 7 (the string bispectrum amplitude 
reduced by a factor of 8). For the same string simulation, the 
shape correlations for different Gyu/c 2 were 99.4% or above and 
the amplitude scaled as expected with (G/i/c 2 ) 3 to within 2%. 
The string bispectrum shown in Fig.[TT]is well converged with 
random errors from the averaging procedure small relative to the 
actual signal. We conclude that, assuming the physics and nu- 
merical accuracy of the string simulations that are available, we 
have extracted a string bispectrum of sufficient accuracy for the 
present non-Gaussian analysis. 

The overall three-dimensional reconstruction of the string 
bispectrum shape is shown in Fig. [12] normalized in the usual 
way to approximately illustrate the signal-to-noise expected 
(that is, removing an overall /" 4 scaling, by dividing by the con- 
stant Sachs- Wolfe bispectrum shape). The first point to note is 
that the bispectrum is negative, reflecting the underlying string 
velocity correlations and curvature correlations that have cre- 
ated it; in the expanding universe, curved strings preferentially 
collapse, creating a negative temperature fluctuation towards the 
centre and a positive signal outside. In the overall spectrum, the 
n=0 mode is dominant, but it is modulated by other modes pro- 
viding further interesting structure that could be described as 
broad arms extending along each axis (see Fig. [TT]); although 
somewhat "squeezed", the correlation with the local model is 
low. The string simulation power spectrum shown in Fig. [5] can 
be understood to be quantitatively modulating the string bispec- 
trum away from the axes, with the signal slowly decaying be- 
yond (say) h,h> 500 in the I3 direction. 

The CMB bispectrum and trispectrum induced by the post- 
recombination gravitati onal effects of cosmic strings have been 
estimated anal ytically ( [Hindmarsfi et al.||2QQ9] |2010[ |Regan & 
Shellard 2010| ). With simplifying assumptions, these predicted 
that the constant mode would be dominant with a broad central 
"equilateral" peak, but not the substructure observed in Fig. [12] 
In terms of missing physics in this post-recombination string bis- 
pectrum, we expect the recombination signal to lie in the range 
I = 200 - 1000 (shown in the full NAMBU power spectrum in 
Fig. [3]) and to significantly enhance the the overall am plitude of 
the bispectrum (see also Landriau & Shellard 2011, where re- 
combination physics is included). 

The correlation of the post-recombination string bispectrum 
with standard primordial shapes is small, because it does not 
contain an oscillatory component from the transfer functions. 
The local, equilateral and orthogonal bispectrum models corre- 
late with strings at about 6%, 11% and 12%, respectively. The 
late-time CMB ISW-lensing bispectrum b 1 ™^ is also mainly 
negative, but it is much more squeezed/flattened and correlates 
at only about 1 1 % with the present string bispectrum, so the pre- 
dicted ISW signal should provide a small positive bias for string 
detection of about 0.44 cr (see below). Diffuse point-source con- 
tamination is of greater concern because at £ max = 2000 this 
has a 40% anti-correlation with strings (for the simple Poisson- 
distributed point source template with b^^ = constant). This 
close relationship with point sources requires a joint analysis. 
Other foreground contamination must also be considered, as we 
shall discuss, and for this we rely on realistic simulated fore- 
ground residual s provided by the Planck co mponent separation 
pipelines (see Planck Co llaboration XII 2013]). 




1500 2000 



Fig. 12. Modal reconstruction of the post-recombination string 
bispectrum Eq. (14) extracted from Planck resolution map sim- 
ulations. This is a 3D view of the allowed tetrahedral set of mul- 
tipoles (€\,€2, £3) showing isosurfaces of the bispectrum density 
with darker blue for more negative values (it is normalized rela- 
tive to the constant SW bispectrum). 



4.2.3. Planck string bispectrum results 

Using the modal bispectrum estimator, we have searched for the 
string bispectrum in the Planck CMB maps obtained using the 
foreground- separation techniques SMICA, NILC and SEVEM. We 
note that the modal estimator has passed through the full vali- 
dation suite o f NG tests described in the IPlanck Collaboration! 
XXIV ( 2013| ), where further details about the analysis can be 
found. In summary, we have used the standard U73 mask, which 
includes a Galactic cut and a conservative point source mask, 
together with "inpainting" as in the /nl analysis (essentially 
apodizing the mask). Together with the foreground separated 
maps, realistic noise simulations were were used to determine 
the estimator's linear correction term and to determine the bis- 
pectrum variance, which was very nearly optimal. For calibra- 
tion purposes, we always compare to the string model with ten- 
sion G///c 2 = 1 x 10" 6 , for which can choose to define and nor- 
malise a string bispectrum parameter / NL = 1 at £ max = 2000. 
The standard deviation A/nl = 0.2 obtained from this Planck 
analysis would imply G/i/c 2 = 1 x 10" 6 string detection at 5cr. 
Instead, defining an F NL normalized relative to the local model, 
we could expect to measure F^l = 31.6 + 6.3. The strong scaling 
of the bispectrum amplitude on the string tension oc (G/i/c 2 ) 3 , 
implies a given measurement yields Gji/c 2 = (/nl) 1/3 X 10" 6 . 

The results of the string bispectrum estimation for each of the 
SMICA, NILC and SEVEM maps are shown in Table]?] Given that 
the Planck data exhibit significant detections of both ISW lens- 
ing and residual point source signals, we also quote their mea- 
surements. In the first instance, we have undertaken an indepen- 
dent analysis of each map, which showed no evidence for a cos- 
mic string signal (all estimates were were within lcr). However, 
given the significant measurement of diffuse point sources and 
their strong anti-correlation with the string bispectrum, we have 
also undertaken a joint analysis of the Planck data (which in 
this case is the same as marginalising over the point source sig- 
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Table 4. Modal bispectrum analysis of foreground-separated 
SMICA, NILC and SEVEM maps showing /nl from strings, ISW- 
lensing and diffuse point sources. Three values for / NL are given 
from independent analysis, joint point source/string analysis af- 
ter ISW-lensing subtraction, and joint analysis after both ISW- 
lensing and foreground residual subtraction. Resulting 95% con- 
fidence limits for Gyu/c 2 are also given. 



Method 


Bispectrum 
signal type 


Independent 
analysis /nl 


ISW subtract 
Joint /nl 


ISW/FG res. 
Joint /nl 


SMICA 


Lensing ISW 
Diff.PSxlO 28 
Cosmic strings 
G/u/c 2 (95%) 


0.75 + 0.37 
1.05 + 0.32 
0.19 + 0.20 
8.4 x 10" 7 


1.35 + 0.34 
0.50 + 0.21 
9.7 x 10" 7 


1.40 + 0.34 
0.37 + 0.21 
9.3 x 10" 7 


NILC 


Lensing ISW 
Diff.PSxlO 28 
Cosmic strings 
Gid/c 2 (95%) 


0.91 + 0.36 
1.16 + 0.32 
0.13 + 0.20 
8.1 x 10" 7 


1.44 + 0.34 
0.46 + 0.21 
9.6 x 10" 7 


1.44 + 0.34 
0.23 + 0.21 
8.7 x 10" 7 


SEVEM 


Lensing ISW 
Diff.PSxlO 28 
Cosmic strings 
Gid/c 2 (95%) 


0.6 + 0.36 
1.07 + 0.35 
0.10 + 0.20 
7.9 x 10" 7 


1.33 + 0.38 
0.38 + 0.21 
9.3 x 10" 7 





nal). Before doing so, we subtract the expected ISW lensing sig- 
nal that provides a /nl = 0.09 string bias. The third column 
in Table [4] shows that marginalising point sources enhances the 
string signal up to the 2cr level for all component separation 
methods; essentially the constant mode becomes more strongly 
negative once the measured point sources are removed. 

Other foregrounds, such as dust emission, could poten- 
tially produce spurious string signals if not subtracted prop- 
erly. Foreground contamination has been extensively simulated 
within the Planck pipeline and foreground residual maps are 
provided by each component separation team. We have anal- 
ysed the residual maps provided by both SMICA and NILC to 
seek evidence of correlations with the string bispectrum. The 
SMICA combined-residual map, without point sources and anal- 
ysed with realistic noise, produces a string bias of /nl = 0-23, 
which after ISW subtraction becomes /nl = 0-14 (relative to a 
variance A/nl = 0.20). After both ISW and foreground residual 
subtraction, a joint analysis with point sources yields a SMICA 
string signal /nl = 0.37 ± 0.21 . The apparent string bias in 
NILC from residual foregrounds was even higher /nl = 0.22 (af- 
ter ISW subtraction), meaning a joint analysis obtained /nl = 
0.23 ± 0.21 (see Tableg]). 

We conclude, given our present understanding of point 
sources and foregrounds, that there does not appear to be signif- 
icant evidence for a string bispectrum signal in the Planck nom- 
inal mission maps, so we infer the following post-recombination 
bispectrum constraint on strings (from /nl = 0.30 + 0.21): 



G/i/c 2 < 8.8 x 10 -7 (95% confidence) . 



(17) 



The susceptibility of the string bispectrum to point source and 
other foreground contamination deserves further investigation 
and will require improved characterisation of the diffuse point 
source bispectrum (beyond the simple Poisson model), as well 
as identification of other foreground residuals generating a small 
string bias. 



weaker than that from the power spectrum, it is an independent 
test for strings and the first quantitative string bispectrum limit 
to date. This should be considerably improved in future by in- 
clusion of recombination physics and more precise foreground 
analysis. A comparison with the power spectrum amplitude indi- 
cates the string bispectrum should rise by (2) 3/2 , which, together 
with the full mission data, would see the sensitivity improve by 
a factor of two (allowing constraints around Gji/c 2 < 4 x 10" 7 ). 
We note that the bispectrum is not the optimal non-Gaussian test 
for strings, because the string signal is somewhat suppressed by 
symmetry (the bispectrum cancels for straight strings). This fact 
motivates further study of the trispectrum, for which the Planck 
sensitiv ity is forecast to be AG/u/c 2 « 1 x 10" 7 ( [Fergusson et al 



2010b), as well as joint analysis of polyspectra. 



The string bispectrum constraint Eq. ( 17 ) is a conservative 
upper limit on the string tension Gli/c 2 because we have not in- 
cluded recombination contributions. Although this constraint is 



4.3. Steerable wavelet searches for cosmic strings 

Wavelets offer a powerful signal analysis tool due to their abil- 
ity to localise signal content in scale (cf. frequency) and posi- 
tion simultaneously. Consequently, wavelets are well- suited for 
detecting potential CMB temperature contributions due to cos- 
mic strings, which exhibit spatially localised signatures with dis- 
tinct frequency content. Wavelets defined on the sphere are re- 
quired to analyse full-sky Planck observations (see, for exam- 
ple, |Freeden & Windheuse rjl 997 [ [Wiaux et al. 12005 [|Sanz et al. 
20061 IMcEwen et al.|2006t |Starck et al.|2006[ |Marinucci et al. 
200^ [Wraux et al.|2008| ) 

We perform an a nalysis using the st eerable wavelets on the 
sphere constructed by Wiaux et al. ( 2005 ). Here we exploit steer- 
ability to dynamically adapt the orientations analysed to the un- 
derlying data, performing frequentist hypothesis testing. We ap- 
ply the first (1GD) and second (2GD) Gaussian derivative steer- 
able wavelets, defined on the sphere through a stereographic 
projection, in order to search for cosmic strings in the Planck 
data. A steerable wavelet is a directional filter whose rotation by 
X e [0, 2n) about itself can be expressed in terms of a finite lin- 
ear combination of non-rotated basis filters. Thus, the analysis 
of a signal with a given steerable wavelet *F naturally identifies a 
set of wavelet coefficients, Wy(a)o,x,R), which describe the lo- 
cal features of the signal at each position ojq on the sphere, for 
each orientation x and for each physical scale R. Several local 
morphological properties can be defined in terms of the wavelet 
coefficients ( [Wiaux et al.|2008| ), including the signed-intensity, 

I((Oo,R) = Ww(a)o,Xo,R) • (18) 

This quantity represents the value of the wavelet coefficient at 
the local orientation xo R) that maximizes the absolute value 
of the wavelet coefficient itself. 

The presence of a cosmic string signal in the CMB is ex- 
pected to leave a non-Gaussian signature that induces a modifi- 
cation in the distribution of I(ojq,R) with respect to the lensed 
Gaussian case. We calibrated the dependence of these signatures 
on the string tension using four simulations of the cosmic string 
contribution ( Ringeval & Bouchet [2012| ) combined with a large 
set of lensed Gaussian CMB realizations, along with a realistic 
description of the Planck instrumental properties (refer to Planck 
|Collaboration"XlIl ( |20T3] )). 

A wide range of string tension values were explored, 
Gyu/c 2 e [2.0 x 10" 7 , 1.0 x 10" 6 ], considering several wavelet 
scales, R = [4.0,4.5,5.0,6.0,8.0, 10.0] arcmin. We choose the 
wavelet scale range as a trade off between the signal-to-noise 
ratio of the string contribution and the small scale foreground 
contamination. In fact, the wavelet for the smallest scale con- 
sidered in this analysis peaks at £ = 1300, while extending at 
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Fig. 13. Deviation of the kurtosis of the signed-intensity as a 
function of Gyu/c 2 , normalized to the standard deviation of CMB 
and noise simulations: AK/cr = (K(R,G/i/c 2 ) - K(R,G/i/c 2 = 
0))/cr. The left panel shows results for the 1GD wavelet and 
the right shows the 2GD wavelet. Each curve corresponds to 
a wavelet scale, R (arcmin), included in the analysis. The final 
sensitivity of the method is determined by combining the two 
wavelets and all the scales. 



nent separation pipelines ( Planck Coll aboration XII 2013 ). We 
found that at the smallest wavelet scale considered in this anal- 
ysis, R = 4 arcmin, these residual foregrounds induce a shift in 
the kurtosis, AK/cr = 0.03, i.e., a bias that is negligible for the 
present analysis. This shift increases to more than AK/cr = 0.3 
when extending the analysis to wavelet scales as small as R = 2.5 
arcmin. However, the shift induced by unresolved point sources 
increases the kurtosis of the signed intensity, so it could not rec- 
oncile the tension between the model and the data pointed out 
above. Validation tests performed on the String Challenge data 
set demonstrated that the sensitivity of the estimator can improve 
to AGyu/c 2 « 4 x 10" 7 for small values of G/d/c 2 when wavelet 
small scales are included. For this reason, we believe it could 
be a powerful tool in constraining strings, even though further 
investigation is needed to reliably modeling the string signal, 
the foreground contribution and the noise properties, in order 
to fully exploit the sensitivity of the Planck data. 

Finally, we note that we have also endeavoured to study the 
simulated string maps using spherical wavelets, making an ex- 
tension of previous work ( Wia ux et al.||2008[ |Hammond et al. 
2009) to compute the Bayesian posterior distribution of the 
string tension. Both the spherical and steerable wavelet meth- 
ods offer good prospects for improved non-Gaussian string con- 
straints from the Planck full mission data. 



higher multipoles with a broad distribution. We use maps at an 
HEALPix resolution of N^q = 2048, including multipoles till 
^max = 2500. We analyse the simulations with the same U73 
mask on the Planck CMB map (refer to |Planck Collaboration 
|XII|2013| ), which masks both diffuse and compact foregrounds, 
leaving 73% of the sky remaining for further analysis (refer to 
discussion in Sect. |4. 2 3] ). 

The string non-Gaussian signatures are characterized in 
terms of the kurtosis of the signed-intensity I(coq,R) in Eq. ( [T8] ) 
at the different scales R and for both the 1GD and 2GD wavelets. 
The averaged results from the non-Gaussian simulations were 
used to model the distribution of the kurtosis as functions of 
Gji/c 2 , i.e., K(R,G/i/c 2 ). Other statistics, such as the skewness 
and the Higher-Criticism, have also been explored. We found 
that the kurtosis sensitivity to the string tension is higher than 
the alternative measures. In Fig. [13| we show the difference be- 
tween the average kurtosis at several Gji/c 2 values and the av- 
erage kurtosis for G/u/c 2 = 0, normalized to the standard devi- 
ation of the simulations. On the given range of scales, the 2GD 
wavelet appears to be more sensitive to the string signal. The 
final sensitivity of the method in recovering the string tension 
from simulated data is established through a goodness-of-fit test, 
performed jointly on the two wavelets for all the scales, and tak- 
ing into account the correlations by means of a covariance ma- 
trix estimated from CMB and noise simulations. From the distri- 
bution of the Gji/c 2 values recovered from simulations the null 
hypothesis can be rejected at 95% CL for Gjd/c 2 > 7 x 10" 7 . 

The procedure described above was applied to the SEVEM 
CMB cleaned map, observing no evidence of a string signal at 
any of the scales studied. This would imply a wavelet constraint 
on the string tension at 95% CL, G/u/c 2 < 7 x 10" 7 , however, 
we note that the signed intensity kurtosis of the data is smaller 
than that of the model even for small values of Gji/c 2 , and it 
is not fully compatible with it. A similar behavior is observed 
for the skewness of I(coo,R). We evaluated the impact on the 
estimator of unresolved point sources, both from radio galaxies 
and sub-millimetre star-forming galaxies, using simulations of 
these astrophysical components processed through the compo- 



4.4. Real space tests for cosmic strings 
AAA. Minkowski functionals method 

Minkowski Functionals (MFs) describe morphological proper- 
ties of the CMB field, and can be used as generic estimators 
of non-Gaussianities. They have long been considered to con- 
strain cosmic string physics, for exam ple on gradient tempera- 
ture maps (see e.g., |Gott et aL| [T990). Indeed they have sensi- 
tivity to non-Gaussianity sourced by strings at all orders (i.e., 
including the kurtosis or trispectrum) and they could prove to be 
a powerful tool to constrain topological defects in general. For 
the sake of brevity and conciseness, precise definitions of MFs 
and analytic f ormulations are presented in |Planck Collaboration 
XXIII ( 2013 ) and, here, we only review how MFs can be used to 
constrain th e string energy density Gji/c 2 , following the method 



discussed in Ducout et al. 



( 2012) , for the local model f\ 



local 
NL 



We measure the four normalize^ functionals Vk (k = 0, 3) 
(respectively Area, Perimeter, Genus and Af c i uster ), computed on 
n± = 26 thresholds v, between v m i n = -3.5 and v max = +3.5 
in units of the standard deviation of the map. The curves are 
combined into one vector y (of size n = 104). 

The principle is then to compare measurement of MFs on 
Planck data y to measurements on NG simulations including 
a string tension level G/d/c 2 , that is, we determine y(Gju/c 2 ) 
numerically. These simulations reproduce the same systemat- 
ics (noise, beam) and the same processing (filter, masks) as the 
data. Assuming that MFs are multi-variate Gaussians, we obtain 
a posterior distribution for Gfi/c 2 through a^ 2 test : 



P(Gjd/c 2 \y) ocexp 



X 2 (y,Gfi/c 2 ) 



(19) 



2 Raw Minkowski functionals V* depend on the Gaussian part of 
fields through a normalization factor A k , that is a function only of the 
power spectrum shape. We theref ore normalize functionals v k = Vk/A k 
to focus on non-Gaussianity, see |Planck~C ollaboration XXIII ( |2013| > 
and references therein. 
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with 

X \% Gfi/c 2 ) = [y- y(G/i/c 2 )] T C" 1 [y - y(G^i/c 2 )] . (20) 

The covariance matrix C is computed from 10 4 Gaussian simu- 
lation^Jbecause, given the existing stringent constraints on cos- 
mic strings, this should be accurate without biasing results. The 
cosmic string MF curve y(G/i/c 2 ) is calibrated on 10 3 realistic 
lensed Planck simulations, to which we have added a string com- 
ponent at a specified level. These simulations take into account 
the asymmetry of beams and the component separation process 
(FFP6 simulations, see Planck Collaborat ion ES| ( |2013| ) for a 
detailed description). For the string component, we h ad at our 
disposal only two high resolution string simulations (Ringeval 
|& Bouchet|2012| ), so our model is the averaged curve obtained 
from this combination of Planck and string simulations. 

Due to the nonlinear dependence of MFs on G/u/c 2 and the 
small number of string simulations, the posterior distribution is 
quite complex and noisy. For this reason, we evaluated the pos- 
terior at n NL = 51 values of Gyu/c 2 , between and 10 x 10" 7 , to 
obtain our Planck estimate for Gji/c 2 . This estimate is stable and 
has been validated in realistic conditions with the Planck String 
Challenges described above, and for which we found consistent 
results with the underlying (unknown) Gji/c 2 . 

4.4.2. Minkowski functionals results 

For the constraint on Gyu/c 2 , we analysed the foregrounds sepa- 
rated SMICA map at Af side = 2048 and £ max = 2000, using the 
U73 mask (/ S k y = 73% of the sky is unmasked). The small- 
est point sources holes were inpainted. We applied two specific 
Wiener filters to the map, designed to enhance the information 
from the map itself (Wm) and from the gradients of the map 



14 



(Wbi = 1)Wm). The filters are shown in Fig. 

Additionally, we estimated the average impact of some resid- 
ual foregrounds and secondaries (FG) on Gju/c 2 , using the linear 
properties of MFs and foregrounds models processed through 
the Planck simulation pipeline (FFP6 simulations, see Planck 
Collaboration ES| ( |2013| )). Uncorrected (Poissonian) unresolved 
point sources (PS), Cosmic Infrared Background (CIB) and 



3 The Gaussian simulations endeavour to incorporate realistic noise 
from the Planck data, but only the effective isotropic beam of the com- 
ponent separation method. 




3000 



Fig. 14. The two Wiener filters, Wm and Wm used to constrain 
Gji/c 2 with Minkowski Functionals. 




SMICA, Wu + Wm 



4 6 

Gfi/c 2 (xlO" 7 ) 



8 



10 



Fig. 15. Posterior distribution of the parameter G/u/c 2 obtained 
with Minkowski functionals. This estimate takes into account the 
lensing of the data, but not the effects of foreground residuals. 



Sunyaev-Zeldovich cluster] (SZ) signals can be introduced as 
a simple additive bias A_y p ^"- on MF curves following: 



: /Gsubtracted + A -PS + A -CIB + Ay SZ 



(21) 



These biases are obtained as an average from 100 simulations, 
however, these do not comprehensively cover all the different 
component contributions in the actual Planck data. 

We eventually obtain the posterior distribution of Gfi/c 2 , and 
we integrate it to report confidence intervals. Results are summa- 
rized in Table|5j for raw data (lensing subtracted) and foreground 
subtracted data (PS, CIB and SZ subtracted). The discrepancy 
between the two filters can be explained because the derivative 
filter Wm scans smaller scales than Wm so it is more easily bi- 
ased by foreground residuals. Given the remaining foreground 
uncertainties, we take the most conservative MF constraint for 
the cosmic string contribution to the Planck data to be 



Gji/c 2 < 7.8 x 10 



-7 



at 95% C.L. 



The corresponding posterior is presented in Fig. [15 

Some caveats need to be mentioned that may influence these 
results. First, for the MF method itself, an important limitation 
is the small number of string simulations used to calibrate the 
estimator. The estimator appears to be mostly sensitive to low- 
redshift strings (infinite strings, with redshifts between and 
30), and this is affected by cosmic variance. As low-redshift 
string simulations are much faster to produce than complete 
simulations back to recombination, it should be possible to im- 
prove the robustness of the constraint using these relatively soon. 
Secondly, the impacts of the different point-source foreground 
components (here, PS, CIB and SZ) have been evaluated by av- 
eraging over 100 Planck simulation maps for which the mod- 
elling is only partial. The precise contributions of these differ- 
ent components needs to be investigated in more detail for the 
Planck data. Fortunately, using the linearity of MFs for these 
contributions it will be possible to jointly estimate these as their 
characterisation improves in future studies. Finally, the impact of 
Galactic residuals should also be assessed in further detail, espe- 
cially for the filter Wm that we have observed to be less robust 
against residuals than the Wm filter. 

With advances in studying these experimental effects there 
are good prospects for the full mission data, the sensitivity 



The SZ signal does not include the SZ-lensing NG contribution. 
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Table 5. MFs constraints obtained on Gyu/c 2 , at the 95% C.L. 
These results are obtained on the SMICA map with the U73 mask 
(/sky = 0.73). The "Raw map" result includes only the lensing 
contribution to the data, while the "Foreground subtracted map" 
includes the lensing, Poissonian point sources, CIB and SZ clus- 
ters contributions. 



Gfi/c 2 


W M 


W M 


+ W m 


Raw map 


<6.8xl0" 7 


<7.l 


I xlO" 7 


FG subtracted map 


<6.0xl0" 7 


<3.( 


SxlO" 7 



of the MFs estimator should improve substantially, with sim- 
ulations forecasting possible MF cosmic string constraints of 
Gfi/c 2 < 3 x 1(T 7 at the 95% C.L. We note that further real 
space analysis of string map simulations has been undertaken 
with scaling indices of the pixel temperature distribution (see, 
e.g., |Rath et al.|201 1) . Extensions calculating a set of anisotropic 
scaling indices along predefined directions appear to offer good 
prospects fro string detection. 



5. Conclusions 

We have reviewed the signatures induced by cosmic strings in 
the CMB and searched for these in the Planck data, resulting in 
new more stringent constraints on the dimensionless string ten- 
sion parameter G/u/c 2 . A pre-requisite for accurate constraints 
on cosmic strings is a quantitative understanding of both cosmo- 
logical string network evolution and the effects they induce in 
the CMB. These are computationally demanding problems but 
progress has been made recently on several fronts: First, high 
resolution simulations of Nambu-Goto strings have yielded ro- 
bust results for the scale-invariant properties of string networks 
on large scales, while there has been increasing convergence 
about small-scale structure and loops (for which the CMB pre- 
dictions are less sensitive). Secondly, post-recombination grav- 
itational effects of strings have been incorporated into full- 
sky Planck resolution CMB temperature maps that are impor- 
tant for validating non-Gaussian search methods. Finally, fast 
Boltzmann pipelines to calculate CMB power spectra induced by 
causal sources have been developed and tested at high resolution, 
whether for field theory simulations of strings or textures or for 
models of Nambu-Goto strings. Three-dimensional field theory 
simulations of vortex- strings at sufficient resolution should, in 
principle, converge towards the one-dimensional Nambu-Goto 
string simulations, but testing this is not numerically feasible 
at present. For this reason, we believe it is prudent to also in- 
clude constraints on field theory strings (labelled GjUahX thus 
encompassing cosmic string models for which radiative effects 
are important at late times (such as global strings). We believe 
this brackets the important theoretical uncertainties that remain, 
that is, we have used the best available information to constrain 
both Nambu-Goto strings (NAMBU) and field theory strings 
(AH). This work has also obtained more stringent constraints 
on semilocal strings and global textures. 

5.1. Cosmic string constraints and the CMB power spectrum 

Accurate forecasts for the CMB power spectrum induced by cos- 
mic strings are more difficult to compute than their equivalent 



for simple adiabatic inflationary scenarios. It requires knowledge 
of the source, quantified by the unequal-time correlator (UETC) 
of the defect stress-energy tensor, from well before recombina- 
tion to the present day, which is not computationally feasible. 
Fortunately, we can exploit scale-invariant string evolution to 
extrapolate the results of simulations with substantially smaller 
dynamic range. We use two methods to obtain predictions for 
the UETCs. An unconnected segment model (USM) is used to 
model the properties of an evolving string network, determin- 
ing its density from an analytic one- scale evolution model, and 
the sources are coupled to the line-of- sight Boltzmann solver 
CMB ACT. A second independent pipeline measures the UETCs 
directly from string simulations in the abelian-Higgs field the- 
ory, passing these to a modified form of the CMBEASY Boltzmann 
code. The resulting Nambu-Goto and abelian-Higgs string CMB 
power spectra are illustrated in Fig. [3] Free parameters in the 
USM model can be chosen to phenomenologically match the 
field theory UETCs (denoted the AH-mimic model) and the 
comparison is also shown in Fig. [3] validating the two indepen- 
dent pipelines. 

To compute constraints on cosmic string scenarios we have 
added the angular power spectrum to that for a simple adia- 
batic model, assuming that they are uncorrelated, with the frac- 
tion of the spectrum contributed by cosmic strings being / 10 at 
£ = 10. This has been added to the standard 6 parameter fit using 
C0SM0MC with flat priors. For the USM models we have obtained 
the constraint for the Nambu-Goto string model 

G/i/c 2 < 1.5 x 1(T 7 , /io < 0.015 , (22) 

while for the abelian-Higgs field theory model we find, 

G/i AU /c 2 < 3.2 x 10" 7 , /io < 0.028 . (23) 

The marginalized likelihoods for f\o and in the fio-fl^h 2 plane 
were presented in Fig. [10 With Planck nominal mission data 
these limits are already about a factor of two more stringent 
than the comparable WMAP 7 -year string constraints and these 
Planck limits improve further with the inclusion of high-/' data. 

5.2. Non-Gaussian searches for cosmic strings 

Complementary searches for non-Gaussian signatures from cos- 
mic strings were performed and we have reported constraints 
from the string bispectrum, steerable wavelets and Minkowski 
functionals. These methods participated in the Planck String 
Challenges and have undergone non-Gaussian validation tests. 

The post-recombination string bispectrum has been recon- 
structed and calibrated from string-induced CMB maps using a 
modal estimator. String Challenge analysis with Planck-realistic 
noise simulations and mask indicated a nominal mission sen- 
sitivity of AGfi/c 2 « 5.8 x 10" 7 . Analysis of SMICA, NILC and 
SEVEM foreground- separated maps has yielded /nl = 0-37 ± 0.21 
for the string bispectrum shape, which translates into a bispec- 
trum constraint on the string tension, Gfi/c 2 < 8.8 x 10" 7 (95 % 
CL). Steerable wavelet methods have been calibrated on string 
simulation maps added to Gaussian CMB maps with realistic 
noise and masking, showing a sensitivity of up to AGjd/c 2 « 
4 x 10" 7 . The string signal was shown to have greater impact 
on the kurtosis of the signed-intensity than on its skewness, 
and no evidence of a string signal was found in the Planck 
data. Minkowski functionals have been applied to string sim- 
ulation maps in a Planck-realistic context, computing the four 
functionals — area, perimeter, genus and Af c i uster — after applica- 
tion of Weiner filters. Using these distributions, a Bayesian es- 
timator has been constructed to constrain the string tension. 
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Analysis of the SMICA foreground-cleaned maps yielded a MF 
constraint of G/j/c 2 < 7.8 x 10" 7 (95 % CL). 

Non-Gaussian searches for strings are complementary to the 
power spectrum analysis and yield constraints as low as G/u/c 2 < 
7.8 x 10" 7 , though we note the potential impact of foreground 
residuals in limiting current precision. These are conservative 
upper bounds because they only include post-recombination 
string contributions, unlike the string power spectrum analy- 
sis. Having such a broad suite of tools, ranging from multipole 
space, through wavelets, to real space detection methods, allows 
cross-validation and reinforces the conclusion that there is at 
present no evidence for cosmic strings in the Planck nominal 
mission data. 
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